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A dilute homogeneous 3D Fermi gas in the ground state is considered for the case of a repulsive 
pairwise interaction. The low-density (dilution) expansions for the kinetic and interaction energies of 
the system in question are calculated up to the third order in the dilution parameter. Similar to the 
recent results for a Bose gas, the calculated quantities turn out to depend on a pairwise interaction 
through the two characteristic lengths: the former, a, is the well-known s-wave scattering length, 
and the latter, 6, is related to a by 6 = a — m{da/dm), where m stands for the fermion mass. To 
take control of the results, calculations are fulfilled in two independent ways. The first involves 
the Hellmann-Feynman theorem, taken in conjunction with a helpful variational theorem for the 
scattering length. This way is used to derive the kinetic and interaction energies from the familiar 
low-density expansion of the total system energy first found by Huang and Yang. The second 
way operates with the in-medium pair wave functions. It allows one to derive the quantities of 
interest "from the scratch", with no use of the total energy. An important result of the present 
investigation is that the pairwise interaction of fermions makes an essential contribution to their 
kinetic energy. Moreover, there is a complicated and interesting interplay of these quantities. 
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I. INTRODUCTION 

Recent experiments with magnetically trapped alkali 
atoms significantly renewed interest in properties of 
quantum gases. As it is known, the initial series of these 
experiments concerned a Bose gas {^'^Rb |1(|, ^^Na 0, 
and Li 3]) and resulted in extensive reconsiderations 
and new investigations in the field of the Bose-Einstein 
condensation. In so doing theoretical and experimental 
observations were made that not only confirmed conclu- 
sions made more than forty years ago but also provided 
a new horizon of the boson physics. In particular, one 
should point out good agreement of the results of solving 
the Gross-Pitaevskii equation derived in 1960s Q with 
experimental data on the density profiles of a trapped 
Bose gas P . The so-called release energy measured in the 
experiments with rubidium was also found to be in good 
agreement with theoretical expectations based on the 
time-dependent Gross-Pitaevskii equation 6]. Among 
new theoretical achievements the exact derivation of the 
Gross-Pitaevskii energy functional can be mentioned 
along with the proof that a Bose gas with repulsive inter- 
action is 100% superfluid in the dilute limit As to the 
experimental innovations, observations of interference of 
two Bose condensates ^3 is a good example (for in- 
terestin g th eoretical details see the papers [111, [13 and 
reviews H). 

The first communications concerning experiments with 
trapped fermionic atoms appeared in the literature about 
three years ago when a temperature near OATf was 
claimed to be reached for a trapped ^Li, where Tf is the 
temperature below which the Fermi statistics is of impor- 
tance. Nowadays the temperatures close to 0.2T/ 16] and 
0-5Tf |0| are reported for the ^Li-vapor. Whereas atoms 
of fermionic '^^K were recently cooled down to 0.3T/ . 
So, the regime of the degenerate Fermi gas is already un- 



der experimental study. In view of this fact, reconsidera- 
tion of the basic aspects of the theory of a dilute uniform 
Fermi gas in the ground state is of importance. 

In the present paper a dilute Fermi gas with repulsive 
pairwise interaction is under consideration. Why the sit- 
uation of a repulsive Fermi gas is of interest whereas the 
s— wave scattering length is negative for ^Li and, 
most hkely, for "^^K [l^ so that a trapped ^Li is con- 
sidered as a good candidate for observation of BCS-like 
transition pjl 123. |2^? The point is that the experi- 
ments can produce (and is now producing) the tempera- 
tures at which the BCS pairing does not occur yet. So, 
at these temperature a Fermi gas with attractive pair- 
wise interaction is close enough in properties (with the 
corrections of 0{{T/Tff), where T/T/ « 0.2) to a re- 
pulsive Fermi gas in the ground state. Of course, with 
one obvious alteration: the positive s— wave scattering 
length should be replaced by a neg ative one in the final 
expressions (see, for example, [3). Thus, the experi- 
ments with magnetically trapped atoms of ^Li and ^^K 
offer exciting possibility of exploring the both superfluid 
and normal states of a Fermi gas. In addition, the most 
recent publications [23, 0| demonstrate that there exists 
interesting possibility of ruling the scattering length of 
^Li which varies in a wide range of values, from negative 
to positive ones, when a magnetic field is imposed on the 
system. 

The particular problem to be investigated here con- 
cerns the kinetic i?kin and interaction i?int energies of 
a uniform dilute 3D Fermi gas in the ground state and 
with a repulsive interparticle potential. This problem 
is connected with a more general question related to all 
the quantum gases. The question is if the pairwise in- 
teraction of quantum particles makes contribution to the 
kinetic energy of a quantum gas or not? It is well-known 
that for a classical imperfect gas the pairwise interaction 
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does not make any contribution to the kinetic energy. 
The usual expectations regarding the kinetic energy of 
dilute quantum gases comes from the pseudopotential 
approach. According to these expectations the kinetic 
energy of a dilute quantum gas is not practically affected 
by the pairwise interaction. It means that taken in the 
leading order of the expansion in the dilution parame- 
ter, the total system energy of a ground-state Bose gas 
coinsides with the interaction one if calculated with the 
pseudopotential (see Refs. IHIlllli). For the Fermi 
case the same approach dictates that the kinetic energy 
does not include terms depending on the pairwise po- 
tential in the leading and next-to-leading orders of the 
dilution expansion (see Ref. 24] and Eqs. H24I) and (|25|1 
below). In other words, in what concerns relation be- 
tween the kinetic and interaction energies, a quantum 
gas is very similar to a classical one from the pseudopo- 
tential viewpoint. However, this result was proved to be 
wrong. An adequate and thorough procedure of calcu- 
lating iSkin and Eint of a cold dilute Bose gas has been 
recently developed in Refs. '2^, "2^. It proves that the 
pairwise interaction in a Bose gas has a strong effect on 
the kinetic energy. Moreover, there are quite real situ- 
ations when the kinetic energy of a uniform dilute Bose 
gas is essentially more than the interaction one! It is 
now necessary to clarify this situation in the Fermi case. 
The more so, that the interaction and kinetic energies 
of imperfect trapped quantum gases are now under ex- 
perimental study |20l |24|. Thus, the aim of the present 
publication is to generalize the procedure developed in 
|25l |2g to the Fermi case. 

The paper is organized as follows. The Section ^1 
presents the kinetic and interaction energy of a ground- 
state repulsive Fermi gas found with the HeUmann- 
Feynman theorem on the basis of an auxiliary variational 
relation given in Refs. |253^J. The Section UTTI is to con- 
sider the derived expressions in various regimes: from a 
weak coupling to a strong one. This is needed to discuss 
the failure of the pseudopotential approach in operating 
with E'kin and E'int- Derivation of i^kin and i^int in Sec- 
tion 2 is simple but rather formal so that some questions 
can remain. This is why Sections IIVI and Ivl give a more 
physically sound way of calculating the kinetic and inter- 
action ener gies . This way invokes a method developed in 
the papers [2^ |2g and dealing with the pair wave func- 
tions, which allows one to go in more detail concerning 
the microscopic features of dilute quantum gases. 



with the pairwise interaction V{r) — 7$(r), where 7 > 
is the coupling constant and ^(r) > stands for the 
interaction kernel (r = Irl ) . Below the particle spin is 
assumed to be s = 1/2 [23 • The ground-state energy of 
the system in question E ~ {0\H\0) obeys the well-known 
relation 



SE = {0\SH\0) 



(2) 



called the Hellmann-Fcynman theorem, SE and SH being 
infinitesimal changes of E and H. An advantage of this 
theorem is that it yields important relations connecting 
the total ground-state energy E with the kinetic E'kin and 
interaction Eint energies. These relations read 



Bm \ ^ 9 m ' 



7 



dm 
dE 



0)=E, 



kin ; 



E, 



(3) 
(4) 



If the dependence of the ground-state energy on the cou- 
pling constant and particle mass were known explicitly, 
one would readily be able to calculate Skin and -Eint by 
means of Eqs. Q and I0J. However, it is not the case as a 
rule, and the dependence is usually given only implicitly. 

In the situation of the repulsive Fermi gas the depen- 
dence of the ground-state energy on 7 and m is indeed 
known only implicitly. According to the familiar result 
of Huang and Yang found with the pseudopoten- 
tial approach but then reproduced within the bound- 
ary collision expansion method |29l | beyond any effective- 
interaction arguments, the energy per fermion e = E/N 
reads 
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which is accurate to the terms of order kpO^. In Eq. (|SJ a 
stands for the s— wave scattering length, kp is the Fermi 
wavenumbcr given by 



kp = (37^2n)l/^ 



(6) 



where n = N/V, and the thermodynamic limit N — » 
00, V — > cxD, n — N/V — i- const is implied. Inserting 
Eq. (O in Eqs. © and Q), one can arrive at 



II. HELLMANN-FEYNMAN THEOREM 

Let us consider the system of N identical fcrmions 
placed in a box with the volume V and ruled by the 
following Hamiltonian: 
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where £kin = Ek^n/N and Sint = E-mt/N. Hence, to derive 
the kinetic and interaction energies from Eq. Q with the 
help of the Hellmann-Fcynman theorem, we should have 
an idea concerning the derivatives of a with respect to 
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the particle mass m and coupling constant 7. As E - 
Ekin + Eint, then from Eqs. (TJl and © it follows that 



m da/ dm — 7 da/dj. 



(9) 



This property of the derivatives becomes clear if we re- 
mind that in the 3D case the s— wave scattering length 
is given by 



(10) 



where f{r) obeys the the two-body Schrodinger equation 
in the center-off-mass system: 



[n^ /m-f)V'^ip{r) + $(r)(^(r) = 0. 



(11) 



The pair wave function ^p{r) represents the zero- momen- 
tum scattering state, and (p{r) — > 1 when r 00. The 
scattering part of the pair wave function given by the 
definition (p{r) = 1 + ip{r) is specified by the following 
asymptotic behavior: 



ip{r) —a/r (r — > 00). 



(12) 



Note once more that the pairwise potential involved in 
Eqs. (Cnil and l^lj is V{r) = 7$(r) but not <I>(r) which 
is the repulsive interaction kernel. As it is seen from 
Eqs. and Ijlll) . the scattering length depends on the 
particle mass and coupling constant through the product 
7717. Hence, to use Eqs. ((TJ and JSjl we should know the 
derivative of a with respect to 7717. 

This derivative can be found with very useful varia- 
tional theorem proved in the papers |25l |26| . After small 
algebra the result of this theorem is rewritten in the fol- 
lowing form: 



Sa = / d^r^^ir^r), 



(13) 



where, remind, (/?(r) is a real quantity. In view of crucial 
importance of this theorem, let us make an explaining 
remark concerning the proof. The key point here is to 
represent Eq. (fllHl as 



d^rip^{r)<i>{r) + — I d^r\V^{r)\\ (14) 



which is realized with the help of Eqs. H12|) and 

V(V'VV') = W V?/; + V'V^V'- So, from Eq. (O one gets 



mda/dm — ^da/dj — a — h, 



(15) 



where the additional characteristic length 6 > is of the 
form 



(16) 



Emphasize that b can not be represented as a function 
of a in principle, and the ratio b/a depends on a particu- 
lar shape of a pairwise potential involved. Now we need 



nothing more to calculate the kinetic and interaction en- 
ergies of the uniform repulsive Fermi gas in the ground 
state. Equations Q and (jS)) taken in conjunction with 
Eq. (|15|l . result in the following expressions: 
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whose sum is, of course, equal to Eq. Q. We again have 
series expansions in kpa but with coefficients depending 
on the ratio b/a. 



III. FROM WEAK TO STRONG COUPLING 

To go in more detail concerning Eqs. (|17|l and H18|l . let 
us consider them in various regimes. We speak about the 
week coupling when the interaction kernel <i>(r) is inte- 
grable and the coupling constant 7^1. The integrable 
kernel with 7^1 and a singular pairwise interaction 
like the hard-sphere potential are related to the strong- 
coupling regime. The expansion parameter kpa involved 
in the expressions mentioned above corresponds to the 
dilution limit fc^ ^ 0. In this situation one is able to op- 
erate with Eq. (jSJ in the both weak- and strong-coupling 
cases. However, for the weak coupling kpa \& small even 
beyond the dilute regime due to a cx 7 ^ 1. This is why 
Eq. (O can be used and rearranged in such a way that 
to derive the weak-coupling expansion for e. 

In the weak-coupling regime the scattering length a is 
given by the Born series: 



a = oo + ai + 
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where Tk = k'^ / {2m) , and <I>(fc) is the Fourier trans- 
form of the interaction kernel ( for more detail see 
Ref. '3^). Inserting Eq. in Eq. 0, one gets the 
following expression: 
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where terms of order 7^ are ignored. Due to Eq. H20|l 
the dependence of Eq. H21|l on the particle mass and cou- 
pling constant is known explicitly. Hence, one can readily 
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employ the Hellmann-Feynman theorem that, taken to- 
gether with Eq. H21|) . yields 
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So, the derived results suggest that the pairwise interac- 
tion influences the both kinetic and interaction energies 
of a Fermi gas. In the weak-coupling regime the major 
part of the 7-dependent contribution to Eq. 1)2111 is re- 
lated to E'int, this part being proportional to 7. While 
the terms of order 7^ appear in both Skin and Eiat ■ This 
conclusion meets usual expectations according to which 
the contribution to the mean energy coming from the 
pairwise potential is mostly the interaction e nerg y for di- 
lute quantum gases (see, for example, Refs. 1^, i§\ and 
the discussion in Introduction of the paper ^Sj). On 
the contrary, beyond the weak-coupling regime the situ- 
ation with i?kin and Ei^t turned out to be rather curious 
and differs significantly from that of the weak-coupling 
case. However, before any detail let us discuss the pseu- 
dopotential predictions for E'kin and Ei^t being the basis 
of usual speculations involving the kinetic and interac- 
tion energies of quantum gases beyond the weak-coupling 
regime. 

At present the customary way of operating with the 
thermodynamics of a dilute cold Fermi gas with repulsive 
pairwise potential is based on the effective-interaction 
procedure: one is able to use either the matrix for- 
mulation like in the Galitskii original paper |3lj | or the 
pseudopotential scheme applied in the classical work of 
Huang and Yang In the dilution limit the t— matrix 
is reduced to t = Airna/m, which yields the momentum- 
independent result inh^a/m for the Fourier transform 
of the effective interaction. This is why one is able not 
to make essential difference between these two effective- 
interaction formulations both referred to as the pseu- 
dopotential approach here. The key point of this ap- 
proach is that to go beyond the weak-coupling regime, 
one should replace the Fourier transform of the pair- 
wise interaction $(fc) by the quantity A-Kh^a/m in all the 
expressions related to the weak-coupling approximation. 
In so doing, some divergent integrals appear due to ig- 
norance of the momentum dependence of the matrix. 
Indeed, substituting t = inh^a/m for $(fc) in Eq. ifHI) . 
one gets a divergent quantity ai oc J d'^k/2Tk that makes 
contribution to the total energy of the system. To de- 
rive the classical result of Huang and Yang, the diver- 
gent term proportional to ai should be removed, which is 



usually fulfilled via a regularization procedure simulating 
the momentum dependence of the i— matrix. In the pseu- 
dopotential scheme of Huang and Yang this corresponds 
to use of the effective interaction {4:nh^a/m)S{r){d/dr)r 
rather than {4:71 fi^ a /m)S{r). From this one can learn that 
to generalize Eqs. ((2U, and to the situation of 
a finite coupling constant, one should replace oq by a 
and remove all the terms depending on ai in the men- 
tioned equations. This yields Eq. |SJl and the following 
pseudopotential predictions for the kinetic and interac- 
tion energies: 



£kin ^ 



10m 



1 - 



2l7r: 



:(ll-21n2)fc|a^ 



-(pb) 



21.2 r 



lOm 



10, 8 

— kpa -\ 

97r 2l7r2 



(11 



-21n2)fc|a2 



(24) 
(25) 



Note that these results can be derived in another way as 
well. For example, the first term in Eq. I|25() can readily 
be reproduced with the pseudopotential in the Hartree- 
Fock approximation (see Ref. [28] and the next section of 
the present paper). From Eqs. and (^31) one could 
conclude that the second term in Eq. Q is related to 
the interaction energy, and, hence, the contribution to 
the mean energy of a dilute cold Fermi gas coming from 
the pairwise potential is mainly the interaction energy. 
However, now we know that actually it is not the case. 
So, one should be careful with the pseudopotential pro- 
cedure which has serious limitations in spite of the cor- 
rect result for the mean energy. Here it is worth noting 
that the pseudopotential scheme preserves some features 
of the weak-coupling regime even being applied in the 
strong-coupling case. This concerns the relation between 
the kinetic and interaction energy for both a dilute Fermi 
ground-state gas and a Bose one ,25, 26] . The same prob- 
lem appears when the pseudopotential is used to calcu- 
late the two-particle Green function in a Bose gas, which 
manifests itself in abnormal short-range boson correla- 
tions Similar troubles can also be expected for the 
two-fermion Green function. 

Now let us consider Eqs. (|17|l and H18|l beyond the 
weak-coupling regime, the ratio b/a being of special inter- 
est. We start with the simplified situation of penetrate- 
able spheres that are specified by the interaction kernel 
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if r > tq. 
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Inserting Eq. H26|) in Eq. (|ll|l . one can find 



r2Asinh(ar)/r if r < vq, 

nr) = i , / -f ^ (27) 

1 — a/r II r > tq, 

where = mj^/h^ ($ > 0) and A is a constant. Equa- 
tion H27() taken together with the usual boundary condi- 
tions at r = To leads to 



a = vq 1 — tanh^aro) / (aro) 



(28) 
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and 



b = ro 



— ^3 tanh(aro ) / (arg ) — csch(aro 



(29) 



where csch(a;) = l/cosh^(a;). One can readily check that 
in the weak-couphng regime, when utq oc 7^/^ — > 0, 
Eqs. if^ and if^ are reduced to 



-aVjj cx 7, 
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and, hence, b <^ a. This means that the next-to-leading 
term in the expansion in kpa given by Eq. (jSJ is mostly 
the interaction energy, as it was mentioned above. On the 
contrary, in the strong-coupling regime, when aro — > 00, 
Eqs. ^ and ^ give 
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(31) 



Hence, b/a 1, and the ground-state energy of a dilute 
Fermi gas with the hard-sphere interaction is exactly ki- 
netic! Note that the same conclusion is valid for a dilute 
Bose gas of the hard spheres [s^ ■ 

Another, a more realistic example concerns a situation 
when the interaction kernel combines a short-range repul- 
sive sector with a long-range attractive one. Here we are 
especially interested in a negative scattering length. It is 
usually considered (see, e.g., 33j) that for alkah atoms 
one can employ the following approximation: 
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+00 if r < vq, 
-C/r'^ if r > To- 



(32) 



The scattering length for the pair interaction kernel (|32|l 
is of the form (see Ref ^J) 

a/r, - r(3/4) J_i/4(a;)/ [2r(5/4) Ji/4(a:)] , (33) 

where x — r1/{2r1), Vc — (rwyC/h^)^^'^, whereas J^ix) 
and r{z) denote the Bessel function and the Euler 
gamma- function. It is known that J,y{x) ~ x'^/[2'^r(l + 
v)] for X — > 0. Therefore, Eq. H33|) reduces to a = tq 
in this limit. In other words, when the attractive sec- 
tor is "switched off", we arrive at the hard-sphere re- 
sult discussed in the previous paragraph of the present 
section. For a; > the scattering length is a de- 
creasing function of x with the complicated pattern of 
behaviour specified by the infinite set of singular points 
{x^ , x^ , x^,...}. These points are the zeros of 
Ji/4{x) so that a -co when x — > — and 

(i) 

a —> +00 when x Xoc + 0. In addition, there is also 
the infinite sequence of the zeros of the scattering length 
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.} being the zeros of J_i/4(x). Note 

< Xq^^\ Keeping in mind this informa- 
tion and Eq. (|33|l . we can explore the ratio b/a for the 
pair interaction kernel (|32|l . Equation (|33l) leads to 

-fda/dj = D^/x J-i/i{x)l {2Ji/i{x)) 

-V2/(^Ji%(x))]. (34) 
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FIG. 1: The ratio b/a versus x = r^/(2ro) for the pairwise 
interaction kernel l|32|l . rc = {rwyC /h^Y^'^ . 

with D = ror(3/4)/[23/2r(5/4)]. Note that to derive 
Eq. (|34|l . the useful formula 

Jv+l{x)J^^{x) + Jy{x)j_{yj^i){x) = -2SIT\{'KV)/{'KX) 

should be applied. Equation H34|l . taken in conjunction 
with Eq. lfT3|) . yields 

b/a = 3/4 + l/[7r\/2Ji/4(a;)J_ 1/4(2;)] . (35) 

As it is seen from Eq. (|35|l . in the limit a; ^ we get 
the hard-sphere result b/a = \ (see Fig. The quan- 
tity b (remind that b is always positive!) is finite at 

X = Xq \ while b —> +00 for x —> Xoo ■ In the latter situ- 
ation b goes to infinity in such a way that b/\a\ — > -l-oo 



though I a 

„(*) 



00 for X 



too. Hence, 



„(') 



and 



xio' are both singular points of b/a. Let us stress that 
the zeros of the scattering length in the case considered 
have nothing to do with the weak-coupling regime for 
which, remind, b/\a\ <^ 1. Operating with the kernel 
(|32l) we are not able to reach the weak-coupling regime 
at all because in this kernel is not bound from above. 
Now let us consider the situation of a negative scattering 
length being of special interest in the experimental con- 
text. The scattering length given by Eq. (|33|l is negative 

provided that Xq*^ < x < xS- As it seen from Fig. 
for any of these intervals the ratio b/a has a maximum 

value [6/a]rn'ax, and it decreases while i increases. In par- 
ticular, [6/a]max ~ —5, whereas [6/a]max ~ —12 and 
[&/a]mL « -18 (see Fig. (|TJ). Hence, Eq. §^ turned 
out to make it possible to get some information about 
b/a even without specifying the range of the relevant val- 
ues of X (in spite of the fact that this range is in principle 
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known). Indeed, according to the mentioned above, the 

ratio h/a does not exceed [6/a]max ~ — 5 if the scattering 
length is negative. This suggests that the contribution of 
the pairwise potential to the kinetic energy is much larger 
than the absolute value of the corresponding contribution 
to the mean energy for a normal-state dilute Fermi gas 
with a negative scattering length at temperatures close 
enough to zero! The interaction energy is negative in this 
case and also much larger, if taken in absolute value, than 
the sum of the a— dependent terms in the Huang- Yang 
result. Note that for alkali atoms [s^l one can expect 
that X ~ 10, which means that h/\a\ ^ 20 (see Fig.m). In 
view of the recent results on a trapped Fermi gas, it 
is also of interest to consider behaviour of 6/a in vicinities 



of the special points and x^' . Varying the magnetic 
field actin g on the system of ^Lz— atoms, the authors of 
the paper |2J were interested in the regime of the Fesh- 
bach resonance, for which a — > oo, and in the situation 
when a — > 0, as well. The both variants, as it follows 
from our consideration, are characterized by b/\a\ ^ 1. 
Note that passage to the limit a ^ oo in Eqs. H17|) and 
(|18ll is not correct because it violates the expansion con- 
dition kp\a\ <^ 1. On the contrary, one can set a = in 
these equation, which leads to the exact result, beyond 
the perturbation theory, 



Skin (a O) 
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(36) 
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Equations (|36|l and H37|) correspond to an unusual and 
extreme situation that, nevertheless, is experimentally 
attainable now (see Ref. jH [21I). The total energy of 
the system is here equal (or practically equal) to that 
of an ideal Fermi gas, while the interaction and kinetic 
energies taken separately have nothing to do with those 
of a gas of noninteracting fermions. The most interesting 
particular case concerns the regime kpb ^ 1, where the 
first term in Eq. is negligible as compared to the 
second one. In this case |ei„t| ~ £kin « n. 

Thus, the examples listed above show that in phys- 
ically relevant situations the correct results for the ki- 
netic and interaction energies of a dilute Fermi gas given 
by Eqs. I|17|l and (|18|l differ significantly (more than by 
order of magnitude!) from the pscudopotential predic- 
tions (|24() and H25() . As it is seen, in a Fermi gas the 
pairwise interaction has a profound effect on the kinetic 
energy contrary to a classical imperfect gas. And this well 
meets the conclusion on an interacting Bose gas derived 
in Refs. (25i. i26l|. A physically sound way of explaining 
this feature of quantum gases is to invoke the formalism 
of the in-medium pair wave functions. This is why below, 
in Sections (|1V|I and 10, the interaction and kinetic en- 
ergies of a Fermi gas are investigated through the prism 
of this formalism. 



IV. INTERACTION ENERGY VIA THE PAIR 
WAVE FUNCTIONS 

The derivation of the kinetic and interaction energies 
of a dilute ground-state Fermi gas given in Section IITI is 
mathematically adequate. However, from the physical 
point of view it has an obvious disadvantage. Namely, 
the microscopic information remains hidden in Eqs. 
H17I) and H18I) due to its implicit usage in Section^ To 
eliminate this shortcoming, £int and £kin are considered 
below with a physically sound approach based on the in- 
medium pair wave functions (PWF). 

For the sake of convenience, let us begin with the inter- 
action energy. It is well-known that all the microscopic 
information concerning the A^— particle system is con- 
tained in the A^— particle density matrix. In the case of 
interest the A— matrix is defined by 



9Nix'i,X2, ■ ■ ■ ,x'^;xi,X2, ■ ■ ■ ,xn) = 
=**(xi,a;2, . . . ,xn)'^{x[,x'2, ■ 



(38) 



where "^{xi,X2^ ■ ■ ■ ,xn) is the ground-state normalized 
wave function, x — {r, a} stands for the space coor- 
dinates r and the spin z— projection a = ±1/2. It is 
also known that actually we does not need to know the 
A^— matrix in detail. In particular, to investigate the to- 
tal system energy together with the kinetic and interac- 
tion ones, we can deal with the 2— matrix defined by 

g2{x[,X2;xi,X2) = jdx^,. . Axn'^* {xi,X2,Xz, ■ ■ ■ ,Xn) 
V 

X ^!{x[,x'2,X3, . . . ,xn), (39) 



where in general 



, = ^ / . ..d 



Let us introduce the eigenfuctions of the 2— matrix 
£,„{xi,X2) given by 

dXidx2 Q2ix'i, x'2; Xi, X2)^i,ixi, X2) = 

=w^^,{x[,x'2), (40) 

where Wi, stands for the v— state eigenvalue. These eigen- 
fuctions are usually called in- medium PWF [s^. The 
2-matrix can be expressed in terms of its eigenf unctions 
and eigenvalues as follows: 

g2{xi,X2;x[,X2) = ^Wi, f*(xi, X2)Cz.(xi, 4), (41) 

V 

where it is implied that 

dxidX2 Q{xi,X2)^„'{xi,X2) = S^^>. 
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From Eq. (jSHJl it follows that 



J dxidx2 g2{xi,X2;Xi,X2) = 1, 



and, hence, 



(42) 



(43) 



which allows one to interprete the eigenvalue as the 
probability of observing a particle pair in the j/— state. 

Now let us remind that the total momentum of 
the system of interest, the total system spin and its 
z— projection are conserved quantities [sH. This means 
that they commute with the iV— particle density ma- 
trix because the latter is permutable with the system 
Hamiltonian. As the total pair momentum HQ, the to- 
tal pair spin S and its Z— component Sz commute with 
the total system momentum, total system spin and its 
Z— projection, correspondingly, they commute with the 
TV— matrix, too. If so, then one can derive that Q, S and 
Sz are permutable with the 2— matrix. This is why we 
can choose the eigenf unctions of the 2— matrix in such a 
way that i' = {^,Q, S,ms}, where ms is an eigen- 
value of Sz and A stands for other quantum numbers. 
Hence, in the homogeneous situation one arrives at (see 
Refs. IH^ and 113) 



(44) 



ii^{xi,X2) = I? i.(r, (71,(72) exp(iQR)/Vy, 



where r = ri — r2 and R = (ri + r2)/2. As the in- 
medium bound pair states like the BCS-pairs are beyond 
the scope of the present publication, here we deal only 
with the scattering states. In other words, only the sector 
of the "dissociated" pair states is taken into considera- 
tion. Hence, — {q,Q, S,mg}, where q stands for the 
relative wave vector. This is why it is convenient to set 
by definition 



(45) 



i?^(r, (71,(72) = fuir,(^i,cr2)/yV 



The pair interaction of interest does not depend on 
the spin variables, which means that (^i/(r, cti, CT2) is ex- 
pressed as 

lfii,{r, (71, (72) = (fi q,Q,s(r) XS,ms{<^l,Cr2), (46) 

where for the singlet states {S ~ 0) one gets 

Xo,o(cti, (^2) = A((7i + (72)sign((7i)/\/2, (47) 
while for the triplet wave functions {S = 1) 

Xl..ms {(^1,(^2) = 

re(-(7i)e(-s2) if ms = -i, 

= < A((7l + (72 )/V2 if ms- 0, (48) 
[e((7i)e(a2) if ms= 1. 



In Eqs. (133 and gHJ 

Now, from Eq. (gT)) it follows that 

Xo, 0(^1,(72) = -Xo,o(ct2,cti). (49) 
Then, the Fermi statistics dictates 

<^q,Q,0(l-) = <y5q,Q,0(-l-) = (P-d^Qfiir), (50) 

and for r — > cxd the wave function (pq,Q,o(r) obeys the 
asymptotic regime 

¥'q,Q,o(r) \/2cos(qr). (51) 
In turn, for triplet states Eq. (|48|) yields 

Xl,ms(0'l)Cr2) = Xl,ms(^2,CTi), (52) 

which leads to 

V'q.Q.iW = -¥'q,Q,i(-r) = -V3-q,Q4(r). (53) 
Now, another boundary regime 

<^q,Q4(r) V2sin(qr). (54) 

is fulfilled when r —^ oo. 

Working in the thermodynamic limit N oo,V 
oo, N/V = n = const, it is more convenient to leave 
the 2— matrix in favour of the so-called pair correlation 
function 

F2{xi,X2;x[,x'2) = {ip'' {xi)-4;'' {x2)ip{x'2)ip{x'i)) , (55) 

where (A) stands for the statistical average of the opera- 
tor A, and t(j^{x), ^l^{x) denote the field Fermi operators. 
The pair correlation function differs from the 2— matrix 
by the normalization factor (see Ref. |35| ) . 

F2{xi,X2;x[,x'2) = N{N ~ l)g2{x[,x'2;xi,X2), (56) 

so that F2 remains finite while g2 approaches zero in 
the thermodynamic limit. Indeed, when V — > 00, N — > 
00, N/V = n ^ const, Eqs. (|41|l and (|56|) . taken in 
conjunction with Eqs. 144|l and (|45|l . yield 

F2(xi,a::2;xi,a::2) = 

[ '^{2^)^ PS.msi^l, QX,Q,s(l-)y'q,Q,s(rO 



S.mc 



X Xs,ms(c^i>cr2)xs.msW,'7^)exp{iQ(R'-R)}, (57) 
where the momentum-distribution function 

pWq,Q)=,lini {N{N-l)wd.Q.s,ms} (58) 

V.,N-^oo ^ ) 
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is finite because 'Wq,Q^s,ms ^ l/V"^ (this follows from 
Eq. when V oo). For ps,ms(Qi Q) one gets the 
relation 



S,ms 



(27r)6 



PS,ms{^, Q) = 



(59) 



resulting from Eqs. H43(l and (jSHl- 

All the necessary formulae are now discussed and dis- 
played, and one can turn to calculations of the interaction 
energy. Using the well-known expression 

Eint^^ J dxidx2V{\ri - r2\)F2{xi,X2;xi,X2) (60) 

and keeping in mind Eq. H57() , one gets the following im- 
portant relation: 



Eint- 



i- / d^r Vir) 
d\d\l 



E 

S,ms 



q,Q) l¥'q,Q,s(r)P, (61) 



(27r)6 

provided the equality 

Xs,ms('^l'^2)XS,ms(0'l,Cr2) = 1 



is taken into account. Equation (|61() directly connects 
the interaction energy per fermion with PWF and, thus, 
with the scattering waves defined by 



and 



V'q^Q^oCr) = '(fq,Q,o(r) - %/2cos(qr) (62) 



■0q,Q,i(r) = <^q,Q,i(i") - \/2sin(qr). (63) 



The scattering waves are immediately related to the 
pairwise-potential contribution to the spatial particle 
correlations. Setting V'q,Q,o(r) — V'q.Q.ili") — Oj oi"; 
other words, ignoring that contribution and taking no- 
tice only of the correlations due to the statistics, one 
arrives at the Hartree-Fock scheme. 

So far we did not invoke any approximation when 
operating with the 2— matrix and pair correlation func- 
tion |38l |. However, taken in the regime of a dilute Fermi 
gas, Eq. (|()T|l can significantly be simplified. Indeed, the 
lower densities, the lower momenta are typical of the sys- 
tem. This means that the pair momentum distibution 
Ps.ms(Q)Q) is getting more localized in a small vicinity 
of the point q = q = when n ^ 0. Consequently, the 
low-momentum approximation can be applied according 
to which for n — s- we get 



d q , ^ \ I / \ 1 9 

g-PS,ms(q:Q) l<y5q,Q,s(r)| ' 



{2n) 



Avsir)\^ Vs,^s, (64) 



where 



and 



(27r)6 



ps.ms (q,Q) 



ips{r) = lim v5q,Q,s(r). 

q.Q^o 



(65) 



(66) 



From Eqs. ^ and ^ it follows that (ps=i{r) = 0. This 
result, taken in conjunction with the low-momentum ap- 
proximation of Eq. lt)4|) . makes it possible to conclude 
that Eq. H61|l reduces for n ^ to 



£int 



2n 



d\V{r)\ipo{r)\^ 



(67) 



The triplet states do not make any contribution to the 
interaction energy in the approximation H64|l . and this 
completely meets the usual expectations. 

Now, to employ Eq. H67() , one should have an idea con- 
cerning 'Po(r) and ?7o,o- As to the limiting wave function 
</?o(r), it can be determined by means of the following 
simple and custom arguments. In the dilution limit the 
pair wave function (^q,Q,s(r) approaches the solution of 
the ordinary two-body Schrodinger equation 

- (?iV"^)VVq,Q,s(r)+T^(r)(^q,Q,s(r) = 

= ih\ym)ipq,Q,sir) (68) 

with the boundary conditions given by Eqs. (|50|l and 
H53(l . Hence, in the limit n —> the quantity (pair) has 
to obey the equation 



(?iV™)VVo(r) + F(0<^o(r) = 0, 



(69) 



where (fioi^^) — > when r — > oo. Comparing Eq. 
with Eq. H69fl . for rt ^ one finds 



(fo{r) = V2ip{r). 



(70) 



To complete calculation of the interaction energy, it only 
remains to find ?7o,o- One can expect that when numbers 
of fermions with positive and negative spin projections 
are the same, the magnitude of ps^ms (Qi Q) appears to be 
independent of the spin variables. In this case Eqs. H59|l 
and give 



VS.,! 



nV4. 



(71) 



Note that Eq. (O can readily be found in a more rigor- 
ous way concerning the relation 



1/2 



d rid r2F2(xi,X2;xi,X2) 



(72) 



where = {iIj^ {xi)4'{xi)). This relation results from 
the definition of the pair correlation function 155f) . To 
derive Eq. (|71|l from Eq. (|72|l . one should employ the 
latter in conjunction with Eqs. H47|l . I|48|l and (|57|1 and, 
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then, take account of n„ = {i/j'^ {xi)tp{xi)) — n/2. Let us 
stress that Eq. H71|l is not general. For example, when all 
the considered fermions have the spin z— projection equal 
to +1/2, one gets rji^i = r? and ?7o,o = ^71,-1 — i]i,o — 0- 
As it is seen, in this case the interaction energy result- 
ing from Eq. (|67|l is exactly amount to zero: one should 
go beyond the approximation defined by Eq. (|64() to get 
an idea about ejnt of such a weekly interacting system. 
Here it is worth remarking that this week interaction is 
an obstacle that can prevent experimentalists from ob- 
serving possible BCS-like pairing of fermions due to an 
extremely low temperature of the BCS-transition. To 
strengthen the interaction effects, it was, in particular, 
suggested 19] to complicate the experimental scheme in 
such a way that fermions with various spin z— projections 
would be trapped. In this case the low-momentum ap- 
proximation yields a finite result for eint- It is deductive 
to go in more detail concerning this situation because 
here another choice of the eigenfunctions of the 2— matrix 
turned out to be convenient rather than that of Eq. 146() . 
The details are in Appendix. 

At last, inserting Eqs. lO and (ITTJ in Eq. and 
making use of Eqs. (|14|l and llt)|) . in the dilution limit 
one can derive eint — 7r7i^n(a — 6)/m, which is nothing 
else but the leading term in Eq. H18|l . Note that to de- 
rive the next-to-leading terms in the expression for the 
interaction energy via PWF, one should construct a more 
elaborated model similar to that of Refs. ^,,26] concern- 
ing a dilute Bose gas. The model like that has to take 
into account in-medium corrections to PWF to go beyond 
the approximation of Eq. Ht)4|l . Though this investigation 
is beyond the scope of the present publication, there are 
some important remarks on the in-medium corrections to 
the eigenfunctions of the 2— matrix in the next section. 

Here it is of interest to check if Eq. (|67|l yields Eq. (|25|l 
when replacing V{r) by the pseudopotential T^(P''^(r). 
The simplest way of escaping divergences while oper- 
ating with the pseudopotential is to adopt V^P^\r) — 
{4:Trti^a/m)S{r) in conjunction with the Hartree-Fock 
scheme. For example, exactly this way was used in 
the classical paper by Pitaevskii when deriving the well- 
known Gross-Pitaevskii equation for the order parame- 
ter of the Bose-Einstein condensation in a dilute Bose 
gas 0. A more elaborated variant, which goes beyond 
the Hartree-Fock framework, requires a more sophisti- 
cated choice of the pseudopotential which, for the par- 
ticular case of the hard-sphere interaction, is of the form 
y(p«)(r) = {AT:h'^a/m)6{r){d/dr)r 2Sj. The aim of this 
variant is to calculate not only the total system energy 
but some additional important characteristics (for in- 
stance, the pair correlation function) which can not be 
properly considered in the former way. Complicating 
the pseudopotential construction allows one to escape a 
double account of some scattering channels (this fact is 
known since the Thesis by Nozieres). This double ac- 
count appears due to the fact that the particle scatter- 
ing makes contribution to the pseudopotential. In our 
case, of course, the simplest choice is enough. Now, re- 



placing V{r) by y(P")(r) = {4Trh'^a/m)S{r) in Eq. ^ 
and setting tp(r) = {(p{r) = 1), for n — > one de- 
rives — {Trh^an/m). It is just the leading term in 
Eq. H25|l . This supports the conclusion that the pseu- 
dopotential is not able to produce correct results for the 
kinetic and interaction energies of dilute quantum gases. 

Here let us make some remarks on the momentum dis- 
tribution of the "dissociated" pairs ps.ms (QiQ)- The cal- 
culational procedure leading to Eq. lt)7|l does not involve 
a detailed knowledge of this distribution. However, it can 
be completely refined. In Ref. |^ it was suggested to de- 
rive ps,ms(Qj Q) via the correlation-weakening principle 
(CWP). According to CWP the pair correlation function 
obeys the following relation: 

F2{xi,X2:,x[,x'2) Fi{xi-x[) Fi{x2;x'2) (73) 

when 

|ri — — > oo, |ri — r'i \ ~ const, |r2 — = const. 

In Eq. (ESIl we set Fi{xi;x'^) = {'^l!^Xl)1p{x[)) . So, the 
pair momentum distribution ps,ms(q, Q), which appears 
in the expansion of F2 in the set of its eigenfunctions, 
can be expressed in terms of the single-particle momen- 
tum distribution n„{k) = (a^(k)ao-(k)), that comes into 
the plane- wave expansion for Fi. In the case of inter- 
est, when the both distribution functions turn out to be 
independent of spin variables, this leads to 

PS,ms(q,Q) =^(IQ/2 + q|)n(|Q/2-q|), (74) 

where, by definition, n{k) — n„{k). Concluding let us 
set, for the sake of demonstration, n{k) = 1 for k < kp, 
n{k) = for k > kp and return to Eq. H64I) . Inserting 
Eq. (|74|l in the right-hand side of Eq. (|64|) and utilizing 
this single-particle momentum distribution of an ideal 
Fermi gas, we arrive at the left-hand side of Eq. (|64|l due 
to /ci? — > when n — s- 0. This example is a good illus- 
tration of the idea of the low-momentum approximation 
introduced by Eq. 

Thus, in Section IIVI it is demonstrated how to calcu- 
late the interaction energy of a dilute Fermi gas from the 
first principles, beyond the formula by Huang and Yang 
taken in conjunction with the Hellmann-Feynman the- 
orem. Though the results of this section make Eq. H18|l 
physically sound and support the conclusion about strong 
influence of the pairwise interaction on the system kinetic 
energy, the nature of this influence is not highlighted yet. 
The detailed discussion of this nature is given in Sec- 
tion |V| 



V. KINETIC ENERGY VIA THE PAIR WAVE 
FUNCTIONS 

Some hints as to how to proceed with the problem of 
the influnce of the pairwise particle interaction on the 
kinetic energy can be found in the Bogoliubov model of 
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a weakly interacting Bo se g as and in the BCS-approach. 
As it shown in Refs. |25ll26j |. there exists some important 
relation which mediates between the pairwise boson in- 
teraction and single-boson momentum distribution nsik) 
in the Bogoliubov model. For the ground-state case this 
relation is the form 



(75) 



where V'-b(^) is the Fourier transform of the scattering 
part of the bosonic PWF corresponding to q = Q = 0, 
and no stands for the density of condensed bosons. When 
the pairwise boson interaction is "switched off" , there is 
no scattering. So, bosonic PWF are the symmetrized 
plane waves and ipBik) = 0. In this case Eq. H75|l has 
the only physical solution nsik) = 0, that corresponds 
to an ideal Bose gas with the zero condensate depletion 
and the zero kinetic energy. On the contrary, "switching 
on" the pairwise interaction leads to ipsik) ^ 0, and we 
arrive at the regime of a nonzero condensate depletion, 
when nsik) ^ and the kinetic energy is not equal to 
zero, as well. 

A similar situation occurs in the BCS-model. There is 
again some corner-stone relation mediating between the 
pairwise interaction and the single-particle momentum 
distribution nscsik). It can be expressed l37| as 

"scs(fc)[l - ?^i3cs(fc)] ^ nlip%csil^)i (^6) 

where ipBCs{k) is the Fourier transform of the internal 
wave function of a condensed bound pair of fermions, Ug 
is the density of this pairs. "Switching off" the pairwise 
attraction leads to disappearnce of the bound pair states: 
"ipBCsik) = 0. In this case there are two branches of the 
solution of Eq. (|75|) : ubcs = 1 and ubcs = 0. Below 
the Fermi momentum the first branch is advantageous 
from the thermodynamic point of view, while the second 
one is of relevance above. So, one gets the regime of an 
ideal Fermi gas with the familiar kinetic energy, often 
called the Fermi energy. When "switching on" the at- 
traction, some significant corrections to the momentum 
distribution of an ideal Fermi gas arise. This corrections 
are dependent of the mutual attraction of fermions and 
make a significant contribution to the kinetic energy ad- 
ditional to the Fermi energy. 

Now, keeping in mind the examples listed above, one 
can suppose that the relation connecting PWF (strictly 
speaking, the scattering waves and bound waves) with 
the single-particle momentum distribution is some gen- 
eral feature of quantum many-body systems. If so, ex- 
actly this relation has to be responsible for the influ- 
ence of the pairwise interaction on the kinetic energy of 
the quantum gases. For a ground-state dilute Fermi gas 
with no pairing effects the relation of interest can be con- 
structed in the form 



n{k)[l-n{k)] = C{k), 



(77) 



where C{k) stands for a functional of the in-medium scat- 
tering waves V'q.Q.s(r). To go in more detail concerning 



the functional C{k), we should try to calculate the kinetic 
energy, starting durectly from n{k) of Eq. H77|l . This 
equation results in 



|[l + VT^T7wl/2 .f ^-^^ 

[[l- ^l-4£{k)\/2 if k>ICF, 

where hICp stands for the Fermi momentum. In the 
present paper a weakly nonideal gas of fermions is under 
investigation ;39], which means that the single- fermion 
momentum distribution approaches the ideal-gas Fermi 
distribution in the dilution limit: C{k) — > and /C^? — > 
kp for n ^ 0. Then, Eq. H78|l can be rewritten for n — s- 
as 

n{k) ~ [1 ~ e{k)] Q{kF -k)+ e{k)e{k ~ kp), (79) 

where the dilution expansions C{k) = £{k){l -I- . . .) and 
ICp = kp{l + . . .) are implied. Taken together with the 
familiar formula 



Eq. (ITU leads for 



to 



3h^kl 2 



d^k 



10m n J (27r)3 



(80) 



(81) 



Note that the characteristic length h given by l|16(l can 
be rewritten as 



27r?i^ 



Sk 
(27r)3 



where ^{k) is the Fourier transform of the scattering 
wave ip(r) (see Eq. (jl2|l '). Keeping this in mind and com- 
paring Eqs. p7|) with H81|l . we can find for n ^ that 



i{k) = (nV4)V^(fc). 



(82) 



Hence, the quantuty C{k) is indeed a functional of 
V'p,q,S'(r) that reduces to the right-hand side of Eq. H82|l 
in the limit n ^ 0. So, our expectations about the re- 
lation mediating between the pairwise interaction and 
single-particle momentum distribution in a dilute Fermi 
gas turn out to be adequate. Let us remark that Eq. ifT^ 
is a good approximation only when calculating the dilu- 
tion expansion for the kinetic energy (strictly speaking, 
the leading and next-to- leading terms). However, to go 
in more detail concerning the single-fermion momentum 
distribution, one should be based on Eq. (|78|l rather than 
on Eq. 179|) . Indeed, Eq. 111|) can be rewritten as 

V'(fc) = j dV<^(r)y(r)exp(-ikr). 

From Eq. (jTUll it follows that / (frLp{r)V{r) = ATih^a/m. 
Therefore, il}{k) oc — 1/fc^ when fc — > 0. Taken in the 
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first two orders of the dilution expansion, the kinetic en- 
ergy is not affected by this singularity. However, n{k) is 
rather sensitive to it. Actually in-medium corrections to 
V'(fc) should be involved to find the adequate behaviour 
of n(fc) when k —^ oo. Indeed, low momenta correspond 
to large particle separations where influence of the sur- 
rounding medium on PWF can be of importance even in 
the dilution regime. 

Discussion about the in-medium corrections to PWF 
can be illustrated as follows. Keeping in mind the re- 
placement of (poi^) by ^/2(p{r) in Eq. H67() . one can sug- 
gest to abandon the approximation C{k) ~ (ri,^/4)?/;^(fc) 
in favour of 



C{k) ~ (nV8)|^o(fc)|' 



(83) 



where ipo{k) is the Fourier transform of ipoi''') — ^oi''') ~ 
\/2. As "00 (fc) V^i^{k) for n ^ 0, the correct dilution 
limit C{k) £{k) results from Eq. H83fl . However, the 
advantage of Eq. H83() as compared with £{k) ~ £{k) is 
that ipoik) is an in- medium PWF and, so, includes in- 
medium corrections at any finite particle density. This is 
why one can expect that Eq. (|83|) . taken together with 
Eq. H78() . gives adequate results for both the kinetic en- 
ergy and the single-fermion momentum distribution. To 
be convinced of this, let us try first to have a guess about 
the in-medium corrections to PWF. Here it is natural to 
start from the expectation that the Fermi sphere is com- 
pletely occupied when n ^ 0. This means (see Eqs. (|77|l 
and (|83|l ') that ipo{k) = ai k < JCp — kp- On the con- 
trary, at large k (small particle separations!) ipo{k) can 
not be significantly affected by surroundings. So, 'ipo{k) 
in this regime is nearly be governed by the ordinary two- 
body Schrodinger equation for the internal wave function 
of the interacting pair at q = 0. Combining these two 
regimes, one can approximately write 

2TkMk) + Qik-kF)J(frifo{r)Vir)exp{-ikr) = 0. (84) 

Stress that the jump in the Fourier transform of ^o{r) at 
k = kp does not imply a jump of i^o{r) itself. Passing to 
the coordinate representation, Eq. H84|l gets the form 



m 



-I' 



Sr'^f^{r')V{r')G{Y-T'), (85) 



where G(r) = /^^^^ rf^fc/(27r)^ exp[ikr]. Equation 
is nothing else but the simplest version of the Bethe- 
Goldstone equation for the internal wave function 
of the interacting pair at Q = q = 0. So, the "veto" 
upon appearance of intermediate scattering states in- 
side the Fermi sphere (see the page 320 in the text- 
book ^35), that allows for passing from the ordinary two- 
body Schrodinger equation to the Bcthe-Goldstone one, 
can be easily explained beyond any intuitive arguments 
with the help of the relation connecting the scattering 



parts of the in-medium PWF with the single-particle mo- 
mentum distribution. 

Now, after the glance at the problem of medium influ- 
ence on PWF, we can try to outline a more elaborated 
way of treating the in-medium pair wave functions. In- 
deed, it turns out that there are everything at our dis- 
posal to derive some two coupled equations which make 
it possible to find ipoik) in conjunction with n{k). The 
first of these equations connecting ipo{k) with n{k) results 
from Eqs. (|78|l and H83|l . As to the second equation, some 
longer but straightforward calculations are needed. The 
point is that Eqs. H67() and H80() enable representation 
of the total system energy as a functional of </?o(?') and 
n{k). Hence, we are able to derive the second equation 
from the minimum condition for this functional. Making 
variation with respect to ipoir) and n{k), one gets 



SE 
'V 



(Pk 
(2^ 



X / d3r(y3o(r)y(r)exp(-ikr) . (86) 



To extract the second equation from Eq. H8t)|) , one should 
keep in mind that the infinitesimal changes Sn{k) and 
SipQlk) are not independent due to the first equation. 
They are related to one another by 

[1 - 2n{k)]Sn{k) = {n'^/8)MkWoik)- (87) 

In addition, a change of n{k) should not affect the to- 
tal number of fermions N = 2 f dPk/{2TT)^n{k). So, the 
equation of interest has the form 



S{E - nN) = 0, 



(88) 



where /i stands for the chemical potential. A combination 
of Eqs. IHEl), 1H3 and (|HHl yields 



Mk) = 



1 - 2n(fc) 
2(m - Tfe) 



d^ripo{r)V{r)e^l>{-ikr). (89) 



Equation H89|) is more complicated but very similar to 
Eq. (|84() . This is especially clear in the light of the fact 
that for n ^ one gets n{k) Q{kp — k). As it is 
seen, when k ^ oo, Eq. (|89|l is reduced to the ordi- 
nary Schrodinger equation for the internal wave function 
of an interacting pair with q = 0. While in the op- 
posite regime, for fc — > 0, Eq. I|89|l acquires completely 
different form with significant contribution of medium- 
dependent terms. So, in-medium corrections are indeed 
of importance for getting an adequate behaviour of n(fc) 
and V'o(^) at small absolute values of the wave vector 
k. It is worth noting that the equation derived by Gal- 
itskii [sJl for the scattering part of his "effective wave 
function of two particles in medium" (see Eq. (11.39) in 
Ref. 01) exactly reduces to Eq. H89|l in the case when 
the center-of-mass and relative wave vectors are equal to 
zero. The author of Ref. (sj introduced the name "effec- 
tive wave function" only due to similarity of the equation 
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for that quantity to the Schrodinger equation for two par- 
ticles in free space. To the best knowledge of the present 
authors, Galitskii did not associate those effective wave 
functions with the eigenfunctions of the 2— matrix. How- 
ever, the derived result strongly suggests that they are 
actually the eigenfunctions of the reduced density matrix 
of the second order. If so, there is a promising possibility 
of using the Galitskii equation in combination with the 
2— matrix, which can produce correct results for e and 
aslo for £kin and eint- Note that in the case of a dilute 
Bose gas the pseudopotential approach does not yield 
correct picture of the short-range boson correlations in 
addition to the failure with the kinetic and interaction 
energies |25| . This is why one can expect the same fault 
with the pseudopotential in the Fermi case. In view of 
this fact, it is worth displaying one more advantage of the 
formalism of PWF. It produces the correct picture of the 
short-range spatial correlations. Indeed, from Eq. (|57|l 
it follows that the pair distribution function g(r\ de- 
fined by g{r) = iMr?)Y.a^,o^-2.{^^^^i\xx,X2), can be 
expressed as 

-^J^^ P5,™.(q,Q) |^q,Q,s(r)r (90) 

When using the low-momentum approximation together 
with Eqs. and ijTlJ, one gets that for r <C l/kp and 
n ^ Eq. (IHOl reads 

g{r) ~ il/2)^\r). (91) 

This result differs from the pair distribution function of 
a dilute Bose gas by the factor 1/2 appearing 

precisely due to the Fermi statistics. Note that this factor 
manifests itself in the total system energy, too. Indeed, 
it is well-known that the leading term in the dilution 
expansion for the total energy of a ground-state uniform 
Bose gas is twice more than the first a— dependent term 
in the corresponding expansion for a Fermi gas. Thus, 
contrary to thepseudopotential calculations for a Bose 
gas (see Ref. [23 ), there is no negative values of g{r) 
at small particle separations when using the formalism 
of the 2— matrix. At last, we remark that it would be 
also of importance to work out some arguments which 
would make it possible to refine the form of £(fc), say, the 
"imperfection" functional. This is of interest in view of 
moving beyond the dilution regime in the investigations 
of the quantum many-body systems. 

VI. CONCLUSION AND COMMENTS 

Concluding, let us highlight the most important points 
of this article. 

First, the low-density expansions for the kinetic and 
interaction energies of a uniform ground-state Fermi gas 
with a repulsive pairwise interaction have been calcu- 
lated up to the third order in the dilution parameter 
kpa. These quantities turn out to depend on the inter- 
particle potential trough the two characteristic lengths 



a and b given by Eqs. (|10|l and H16|) . In the first or- 
ders the 6— dependent terms are cancelled in the sum of 
-Ekin + Emt, and, so, the total energy E, if taken in those 
orders, involves only the scattering length a, as it is well- 
known. However, all the higher orders of the dilution ex- 
pansion for E are expected to emonstrate &— dependence. 

Second, the calculations have been fulfilled in two ways 
with the aim of controlling our conclusions and making 
the present investigation more physically sound. One of 
those ways is based on the Hellmann-Feynman theorem 
which is utilized in conjunction with the useful varia- 
tional theorem for the scattering length when deriving 
i?kin and -Eint from the familiar result by Huang and 
Yang for the total energy of a ground-state dilute uni- 
form Fermi gas. Another way invokes a formalism of the 
in-medium pair wave functions. This variant allows one 
to find -Ekin and i?int starting directly from the first prin- 
ciples. The advantage of the latter way is that it makes 
the underline physics more clear and enables to go in 
more detail concerning the spatial fermion correlations 
and their infiuence on the quantities of interest. 

Third, the ratio of the two characteristic lengths b/a 
has been investigated for the model pairwise interaction 
often used in the context of the alkali-metal atoms. In 
particular, according to the found results, one faces a 
rather interacting gas of fermions in the situation when 
a = 0. The matter is that the sum of the terms coming 
from the pairwise interaction is here exactly cancelled in 
the total energy but, nevertheless, they make significant 
contributions to Ekin and Ei^t , taken separately. This is a 
consequence of presence of the two characteristic lengths: 
when a is amount to zero, then b can take any positive 
value depending on a particular form of a pairwise inter- 
action. Moreover, b can be extremely large in such an 
"uninteracting" gas. Thus, one should be very careful 
with an experimantal analysis of a spatial expansion of a 
trapped gases with a = (see Ref. |23,|21|). Actually in 
such situations the interaction energy can be extremely 
large in absolute value and has the negative sign. 

Forth, important arguments have been presented that 
any many-particle quantum system is characterized by 
some cornerstone relation connecting the single-particle 
momentum distribution with the in-medium pair wave 
functions. Exactly it is responsible for infiuence of 
the pairwise interaction on the kinetic energy of quan- 
tum gases. The form of this relation for a repulsive 
ground-state uniform Fermi gas has been established in 
the leading-order in the dilution parameter kpa. It is 
worth noting that if incorporated in the formalism of the 
2— matrix, the relation under question allows for getting 
Schrodinger-like equations for the in-medium pair wave 
functions. The simplest variant of these equations re- 
duces to the Bethe-Goldstone one. More elaborated ver- 
sion demonstrates important and deep parallels with the 
Galitskii consideration. Thus, one can expect that the 
failure of the pseudopotential approach with i?kin and 
Eint and, also, with the spatial particle correlations oc- 
curs somewhere during its second step when the effective 
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APPENDIX A 

In the appendix the situation of a dilute gas of fermions 
with the single-particle spin ^/2 (£ > 1 is some integer pa- 
rameter) is under investigation. Let particles with some 
two z— projections of single-particle spin be present, for 
example, ^/2 and i/2 — t' , where is integer and £' < £. 
And let the density of the fermions with the projection 
is ni/2, whereas nf/2-£' corresponds to ij^ — l'. This 
situation is of special interest in view of the proposal of 
Ref. 10] where a spin pattern like this was suggested to 
strengthen the pairwise-interaction effects in a magnet- 
ically trapped vapour of ^Li— fermions. In the case un- 
der consideration it turned out to be more convenient to 
compose a set of the eigenfunctions for the 2— matrix in a 
way slightly different with respect to that of Section Hvl 
The difference is that we now adopt v = {q, Q, Sz^ ^Sz\^ 
where Sz stands for eigenvalues of Sz-, the Z— projection 
of the total pair spin, while ^Sz does not correspond to 
some important physical characteristics but simply enu- 
merates degenerate states (for more detail, see formulae 
below). So, here we do not care about the total pair 
spin, in order to simplify all calculations. This choice is 
quite under the general rules of the quantum mechanics 
together with that of Section Hvl (see the discussion con- 
cerning the eigenstates of the 2— matrix in this section). 
Equation 146|l in the situation of interest is of the form 



^^(r, CTi, 0-2) = V^q^qW Xsi(cri, cr2), 



(Al) 
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where the spin variable ct; {i = 1,2) takes the values 
£/2 and £/2 - i' . In Eq. (jJU the condensed notation 
S = 5sz is introduced, and S is taken as a superscript 
to stress the fact that S is here some auxiliary quantity 
rather than an important physical characteristics. The 
eigenstates for Sz are now chosen as follows: 



xriL'>=^(i^>i^-o-io-^')i^>). 



\/2 V'2''2 

\Xsz=e-2e'/ — I 2 ''2 



-n- 



(A2) 



As it is seen, there is no degeneracy for Sz = i and 
Sz — £ — 2t when the only variant Ssz = 1 is involved. 
While for Sz ^ £ — £' two possible eigenstates of Sz are 
present: 5i~gi = 1,2. From Eq. \K2\ it follows that 



and, due to the Fermi statistics, we find 



(A3) 



(A4) 



the wave function 'Pq^Q(r) approaching -\/2sin(qr) when 
r — » 00. On the contrary, for 5 = 2 one gets 



which results in 



(AS) 
(A6) 



where ¥'q^q(i") tends to \/2cos(qr) for r ^ cxi. 

Now, one has everything at his disposal, to express the 
pair correlation function i^2and interaction energy Ei^t in 
terms of PWF. Taken in conjunction with Eqs. ||1T|| . 
and (|45|) . Eq. (|A1I) in the thermodynamic limit yields 

, \ f d\dXl (5), , 



The low-momentum approximation (for more detail, see 
Eq. Ht)4|) ') makes it possible to extremely simplify this 
expression if the dilution limit is of interest. So, for 

11^/2, ni/2-i' one gets 



(2) 



where 



and 



(5) 
^Sz- 



d\(PQl (5), , 



,(^)(r)= lim^^(';)Q(r). 



(AlO) 



(All) 



Only the terms corresponding to 5 = 2 make contribu- 
tion to the interaction energy HA9|I in the dilution limit 
because (p^^\r) is exactly equal to zero as it is seen from 
Eqs. I|A4I) and lAllI) . Following the arguments of Sec- 
tion ^3 (see the discussion just after Eq. one can 
rewrite Eq. (|A9|I in the form 



(A12) 



with (p{r) obeying Eq. pi(l . 

To have an idea about ?7^^^£/, let us turn to Eq. (|72|) . 
Inserting Eq. ljA7p in Eq. H72() . one can find 



E 

Sz.S 



'W) 



pg(q, Q) Ix^°j(ai, a2)|" = n,,n., (A13) 



provided that the necessary normalization condition 



(1/V^) /y c?'^''|</'q*Q(r)P = 1 is taken into account ( 



see 



Eqs. (113 and ijAlJ). From Eqs. I|A2|I and l|A1.3|l it follows 
that 



1 r d\d^ 

2 J 72^ 



P^V(q>Q) + pr-V(q.Q) 



(2) 



»^f/2"-^/2-/!'-(Al4) 



^ \'PqMi^)\^ 1X5^(^1, (72) (A7) Hence, from Eqs. ETO|) and (jXTIjl one can find 



Further, inserting Eq. ljA7|l in Eq. H60|l and keeping in 
mind 

,(-5) 



It is natural to expect that p9\,{q, Q) = pflii^jQ)- 



(A15) 



(■5) 

Ve^t = %/2"£/2-r 



Eixi>i--2)p = i, 



one can arrive at 



Eint/V=^ J d\ V{r 



Ef d\d^ (5) (5) 



(r)l^ (A8) 



Finally, using Eqs. lO, JISIl, l|A12(l and l|A15|l . for 

«£/2, 'ra£/2-£' one gets 

Eint/V ~ [Anh^ia - b)/m] n£/2"£/2-r (A16) 
and, hence, 

E/V ~ {^TTh^a/m) nii2n(,i2~i' , (A17) 

Ey,JV ~ {'^■KhHlm) %/2%/2-£'- (A18) 
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So, as it has been proposed in Ref. [I3 and also seen from 
Eqs. ljA16|l - ljA18|) . trapping of fermions with two differ- 
ent 2— projections of the single-particle spin allows one 
to operate with a system where the effects of the pair- 
wise interaction play a more significant role as compared 
to the situation of extremely weak interacting system of 
one fermion species. This makes it possible to take ad- 



vantage of a large negative triplet scattering length in 
experiments with a degenerate Fermi gas. Indeed, the 
stronger interaction, the higher BCS-transition temper- 
ature And the latter significantly simplifies the ex- 
perimental program of searching for the BCS phase tran- 
sition, especially in the regime \a\ 00. 



